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ABSTRACT 

We present an interpretation of the phenomenological relations between the 
spectral peak, isotropic luminosity and duration of long gamma ray bursts that 
have been discovered by Amati et al., Ghirlanda et al., Firmani et al., and Liang 
& Zhang. In our proposed model, a jet undergoes internal dissipation which pre- 
vents its bulk Lorentz factor from exceeding 1/8 {9 being the jet opening angle) 
until it escapes from the core of its progenitor star, at a radius of order 10 10 cm; 
dissipation may continue at larger radii. The dissipated radiation will be partially 
thermalized, and we identify its thermal peak (Doppler boosted by the outflow) 
with E p k. The radiation comes, in effect, from within the jet photosphere. The 
non-thermal, high energy part of the GRB emission arises from Comptonization 
of this radiation by relativistic electrons and positrons outside the effective pho- 
tosphere. This model can account naturally not only for the surprisingly small 
scatter in the various claimed correlations, but also for the normalization, as well 
as the slopes. It then has further implications for the jet energy, the limiting jet 
Lorentz factor, and the relation of the energy, opening angle and burst duration 
to the mass and radius of the stellar stellar progenitor. The observed relation be- 
tween pulse width and photon frequency can be reproduced by inverse- Compton 
emission at ~ 10 14 cm from the engine, but there are significant constraints on 
the energy distribution and isotropy of the radiating particles. 

Subject headings: gamma rays: bursts; supernovae: general; radiation mecha- 
nisms: non-thermal 
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1. 



Introduction 



The mechanism by which gamma-rays are created in the transient events know as 
gamma-ray bursts (GRBs) has remained surprisingly resistant to interpretation. Two key 
ideas have guided theoretical efforts. The first is that thermalization must be rapid and 
nearly complete close to the engine. A fireball is created, which could reach a temperature 
as high as ~ 1 MeV within a region the size of a neutron star or stellar-mass black hole 
(Paczyriski 1986; Goodman 1986). The second idea is that strongly non-thermal particle 
distributions are expected at large distances from the central engine, e.g. when the relativis- 
ts ejecta collide with the ambient medium (Rees & Meszaros 1992). Continuing dissipation 
within the outflow itself can also be expected, due to the formation of internal shocks (Rees 
& Meszaros 1994), or the release of magnetic energy by reconnection (Thompson 1994; 
Drenkhahn & Spruit 2002) or global current-driven instabilities (Giannios & Spruit 2006; 
Lyutikov k Blandford 2003). 

The prompt emission of a GRB appears at first sight to be highly non-thermal. A 
high-energy cutoff is seen only in a modest fraction of GRB spectra (Pendleton et al. 1997; 
Ryde 2005). Often the low-energy spectrum appears to be harder than would be consistent 
with optically thin synchrotron emission (Crider et al. 1997; Ghirlanda, Ghisellini, & Cel- 
lotti 2003). A discrete thermal bump is usually not apparent in the usual spectral analyses, 
although Ryde (2005) has shown that up to 30% of long GRB can be interpreted as a combi- 
nation of a thermal peak plus a power law spectrum. Intriguing evidence has, furthermore, 
emerged, which indicates that the gamma-ray emission is seeded by radiation with a nearly 
black body spectral distribution. This evidence comes from measurements of the photon 
energy E pk at which the GRB energy spectrum E^dN^/dE^ has a maximum. 

There has been longstanding interest in the distribution of E pk values in GRBs. BATSE 
measured a distribution that is clustered around 200 keV (e.g. Preece et al. 1999). Observa- 
tional selection might explain the high-energy cutoff to this distribution: the GRB spectrum 
is harder below the peak and so a lower flux is measured within the BATSE band from 
bursts with large E pk (Piran & Narayan 1996). Nonetheless, the existence of a substantial 
population of GRBs with E pk higher than ~ 1 MeV remains an open question. 

Our knowledge of the E pk distribution has improved dramatically with the localizations 
provided by BeppoSAX and HETE-II, which have allowed a direct measurement of the 
distances to bursts (see Piran 2005; Meszaros 2006, for detailed references). The value of 
-Epk in the cosmological rest frame of the burst source correlates well the isotropic gamma-ray 
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(Amati et al. 2002; Lamb et al. 2005). Bursts with lower peak energies are therefore 
substantially dimmer. Related correlations have been obtained by correcting E iso for the 
beaming angle, as inferred from the measurement of breaks in the afterglow light curve 
(Ghirlanda et al. 2004; Liang & Zhang 2005), and by measuring the joint correlation of E p k 
with E iso and the burst duration (Firmani et al. 2006). 

The breadth of the intrinsic E pk distribution provides valuable information about the 
emission physics. In the sample of bursts with measured redshifts, E l k generally lies in the 
range 100 keV - 1 MeV (Amati 2006). These values are generally lower than what would 
be expected from thermalization close to the engine. Although a broader distribution of 
-Epk can result from prompt thermalization followed by adiabatic cooling (e.g. Thompson 
1994; Meszaros & Rees 2000, Meszaros et al. 2002), a more plausible explanation is that 
thermalization takes place at greater distances from the engine than was suggested by the 
original fireball models. The relation between E iso and E pk that results from this type of 
distributed heating has been investigated recently (Rees & Meszaros 2005; Thompson 2006; 
Pe'er et al. 2006). The peak energy can remain quite high at large distances from the engine, 
after allowance is made for beaming and bulk relativistic motion. Continuing dissipation in 
the jet is a natural consequence of its interaction with the core and envelope of a Wolf-Rayet 
star. 

Bursts with peak energies less than 50 keV (X-ray flashes) have been detected by Ginga 
and HETE-II (Strohmayer et al. 1998; Sakamoto et al. 2005), but few have measured 
redshifts. Their T 90 distribution is, nonetheless, remarkably similar to that of the harder 
spectrum bursts, which suggests that the low E p ^ is not primarily due to a large cosmological 
redshift. Since the formation of a nearly black body spectrum is possible only out to some 
distance from the engine, one expects that the thermal peak energy must cover only a limited 
range. In fact, as we will see, quasi-thermal models of the prompt GRB emission over ~ 10 
seconds cannot easily be extended below E l pk ~ 50 keV without invoking strong adiabatic 
cooling. Peak energies much higher than 1 MeV also cannot easily be obtained without 
invoking upscattering by relativistic particles. 

It has been noted (Nakar & Piran 2005) that among long bursts in the BATSE catalog 
whose redshift is unknown, there may be up to 25-30% whose low luminosity and hard 
spectra are incompatible with the Amati et al. relation. The implication is that this relation 
may mark out a limiting value of the spectral peak energy at a given apparent luminosity. 
For the quasi-thermal model discussed here, we point out that in fact at low luminosities 
a departure from the Amati et al. relation is expected, which works in the sense of the 
above discrepancy: relativistic jets with large opening angles and relatively low isotropic 
luminosities should produce spectrally harder bursts. 
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In addition to providing seeds for the prompt gamma-ray emission, a thermal photon 
field also can play an important role in accelerating the outflow (Pacyznski 1990; Shemi & 
Piran 1990). Phenomenological scaling relations such as the Amati et al. relation therefore 
provide crucial information about the relative sequence of thermalization and acceleration. If 
thermalization occurs first, at some characteristic radius, then a tight relation between E pk 
and Else follows naturally. Moreover, it has been suggested that the slope of the — E- lso 
relation can be reproduced most easily if the internal heating is driven by non-radial shear 
instabilities near the base of the jet - where its Lorentz factor Y is comparable to the inverse 
of the opening angle 9 - as opposed to radial instabilities in a relativistic flow with Y ^> 1/9 
(Thompson 2006). This raises the possibility that components of the jet which emerge from 
a Wolf-Rayet star with low baryon loading but also a low entropy (e.g. the jet core) will 
attain lower terminal Lorentz factors, and may manifest themselves in the afterglow phase. 

The approach adopted in this paper is two-fold. After an initial outline in §2 of the gist 
of our thermal interpretation of the Amati et al. relation, we begin by working backward 
from the observed correlation between spectral peak energy and isotropic burst energy. In 
§2 we show that thermalization must continue out to a large distance (~ 10 10 cm) from 
the engine in order to reproduce the normalization of the Amati et al. relation. We show 
that thermalization at this radius is reasonable, and derive a condition for the outflow to 
generate enough photons to establish a blackbody spectral distribution (§2.1). In particular, 
the magnetic field must carry at least ~ 10% of the outflow luminosity. The acceleration of 
the outflow is addressed in §2.2, where it is shown that the blackbody photons can effectively 
accelerate the entrained baryons (and magnetic field) up to a limiting Lorentz factor of 100- 
500. In §2.3 we examine the effect on the E pk -Ei SO relation of having thermalization occur 
before or after the acceleration of the jet material. In the second case, E p k has a much 
weaker dependence on E iso than is observed. 

The alternative works much better. The normalization and slope of the Amati et al. 
relation follow immediately from two simple assumptions: first, that the total jet energy is 
regulated to a value close to the net gravitational binding energy of the Wolf-Rayet core and, 
second, that the GRB-emitting component of the jet has r ~ 1/9 at the core boundary. This 
simple model predicts the joint correlation of E pk with E iso and the burst duration tj derived 
from the data by Firmani et al. (2006). The possibility of a range of jet energies is examined 
in §2.4, where it is shown that the Ghirlanda et al. (2004) relation is consistent with an 
approximate scaling Ej ~ 9J 1 . Thermalization can also continue out to the photosphere in 
jets with heavier baryon loadings, and a minimum thermal peak energy of 50 keV (without 
adiabatic cooling) is derived in §2.5. The implications of the E pk -E iso relation of short GRBs 
for our model are outlined in §2.6, where it is shown that a large thermalization radius, 
similar to that of the long GRBs (Rq ~ 10 10 cm) is implied. The pulses of GRBs are 
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generally narrower at higher frequencies (Fenimore et al. 1995), and it is argued in §3.1 
that this is nicely consistent with IC scattering of the seed thermal photons by a power- 
law distribution of electrons (and positrons) at ~ 10 14 cm from the engine. However, this 
explanation for the pulse duration has an important implication: the IC emission by particles 
with relativistic energies must be beamed in the bulk frame. In addition, the distribution 
of particle energies must have a lower cutoff at transrelativistic energies (7 m i n ~ 1). We 
consider the implications for models of particle heating (by turbulence and shocks) in §3.2. 
The paper conclusions with some summary observations in §3.3. 



2. Quasi-Thermal Model of Burst Emission 

We assume here that the flow can be characterized by a mean Lorentz factor T(r) and 
opening angle 9(r) at each radius r. The opening angle can vary with radius if a non-radial 
magnetic field plays a role in collimation and acceleration (e.g. Lynden-Bell 2003; Vlahakis 
& Konigl 2003). An observer viewing the outflow face on sees a luminosity 



2L 

Aso(r) = {6 3 « 1), (2) 



where Lj is the total output of the engine and 9j is the half-opening angle of the outflow. 

A fraction of the outflow energy is assumed to be completely thermalized at a radius 
Ro, where the outflow Lorentz factor is IV In the present context, we will find that the 
temperature T" in the rest frame of the outflow is low enough that thermally-created pairs 
do not contribute significantly to the specific heat. The bulk frame thermal energy density 
is then ~ aT' 4 , where a is the Stefan-Boltzmann constant. An observer at rest with respect 
to the engine sees a temperature 

V2 /1Kr . x 1/4 

(3) 



Tobs " 3° T ° " " bb {-Ro) 



It is reasonable to assume that strong thermalization takes place only at radii comparable 
or less than the radius of the stellar progenitor, R < R*, e.g. that being the limiting 
radius where oblique shocks can cause dissipation. At this radius, one would expect the bulk 
Lorentz factor T (R ) ~ l/2y/39 j where 9j is the jet opening angle at the surface, since above 
this Lorentz factor shear-driven modes in the jet will have no time to grow on a lengthscale 
~ R*9 (in the bulk frame). Thus, using equations ([2]) and ([3]), and assuming T(R) ~ 9J 1 
and Lj = Ej/tj ~ constant (as motivated by Frail et al. 2001, who find Ej ~ const) one 
obtains 

E r cc ^'V"E^ (4) 
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(Thompson 2006). This is essentially the relation discovered by Amati et al. (2002), if we 
consider that Rq ~ i?* ~ constant (e.g., all progenitors of long GRBs are lacking hydrogen 
envelopes) and the dispersion in t l J 4 (10 1//4 for the majority of long GRBs) is smaller than 
that in E]^ (approximately 10 4 / 2 , two decades). 

This simple derivation of the Amati et al. relation involves making some reasonable 
physical assumptions. However, it is more illuminating to work backwards, and use the 
observed relation (pQ) to derive implications for the physics of the model and the progenitor. 
Combining the black body law ([3]) with eq. (CQ) givegj 

^ = 6xl0»4 2 ^g<j; /2 cm. (5) 
i 

One observes that the Amati et al. relation can be reproduced, but only if thermalization 
takes place at a radius Rq that is much larger than the engine radius (in the conventional 
picture where the engine is an accreting black hole or possibly a rapidly-spinning magnetar). 
If the motion of the jet is only mildly relativistic at this point then R is comparable to the 
core radius of a Wolf-Rayet star. 



2.1. Thermalization in a Plasma of Moderate Scattering Optical Depth 

Consider now the consequences of heating an ionized plasma that is optically thick to 
scattering, but optically thin to free-free absorption. Under what circumstances does the 
radiation field relax to a nearly black body distribution? To get a sense of the Thomson 
optical depth t t that is required, we first assume that the radiation spectrum is nearly black 
body, and then examine the self-consistency of this assumption. 

Synchrotron emission by relativistic electrons (and positrons) is generally a copious 
source of soft photons. In spite of this, a minimal magnetic field is required for synchrotron 
emission to supply enough photons to establish a black body gas. We suppose that a fraction 
e sa of the energy density U' of the outflow is transferred to electrons that radiate at the 
synchrotron self-absorption frequency (at a Lorentz factor 7 sa ). A fraction Eb/{sb + %) of 
the energy of these relativistic electrons is released to synchrotron photons (as opposed to 
IC scattering of the ambient radiation field, with a total energy density U' = e^U' from all 
sources). The number of (new) synchrotron photons is 

S synch ~y £B+ J Q^lheB'/meC 1 ' 



throughout, we use the shorthand X = X n x 10", where the quantity X is measured in c.g.s. units. 
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Let us compare this expression with the number density of photons in a black body gas 
(of temperature T bb ) that carries a fraction Ebb of the outflow energy. The characteristic 
self-absorption frequency can be obtained from 



TTsa / B' " ! 



y 5 ^ 

/sa Id 



QED 



(7) 



where TT sa = ^T^'Jcfee)^ is the Thomson depth through the electrons of energy •y sa ,m e c 2 . 
The net output in synchrotron photons at the self-absorption frequency is given by 

4 B' 2 

e sa U' ~ <( 7sa ) x - 7s 2 a0 r T — ct'. (8) 
The optical depth tt S a is therefore regulated to a value 

2 ^sa / n \ 

o7sa T Tsa ~ — ■ (9) 

Combining this expression with eq. ([7j) gives 



V7 / B i \ 



V-B, 



3 QED 



(10) 



The relative numbers of synchrotron and blackbody photons are obtained by substituting 
eq. ([TO]) and the relation {k^T^/ m e c 2 ) [B' j -Bqed) 172 = 1.7 OW^b) 174 in to eq. ©, 



77' 5/7 13/14 / , T , v -3/7 

"■7 synch tsa fcg / ^B^bb \ 

<b(^bb) ~ 4 14 (^+^) ' ^ ^ 

Recall that e sa labels the net energy that is injected into relativistic particles of energy 
^sajrieC 2 . For example, if the energy injected into relativistic electrons is distributed uniformly 
over Lorentz factor and comprises a fraction £ nth of the total, then e sa = e n th/ hi(7max/7min) ~ 
0.l£: n th- The right-hand side of eq. (fTTj) can be larger than unity, but only if 



£bb ~ V £ bb J V m eC 2 



(12) 



We conclude that photon creation in the outflow can be rapid enough to create a black 
body photon gas, but only if the magnetic field carries more than ~ 10 percent of the total 
dissipated energy. 

Soft photons injected into the outflow are redistributed in frequency by multiple Comp- 
ton scatterings. When the electrons and the photons have the same temperature, there is no 
net transfer of energy from electrons to photons, but soft photons will increase their energy, 
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at a rate E 1 /E~ l ~ 4(k-BT^/m e c 2 )aTTi' e c. The efficient redistribution of photons in frequency 
via Compton scattering requires that 



MI 

m e c 2 



TT > 1, (13) 



where tt is the Thomson depth through the warm electrons with a bulk-frame temperature 
T' e . From eq. ffl3|) . one obtains a (conservative) lower bound on the Thomson depth that 
guarantees efficient thermalization of the soft seed photons at an effective temperature T bb , 

TT>( k -^kY. (14) 

We will adopt this value of the optical depth when evaluating the temperature of the radiation 
that emerges from outflows with a relatively high baryon loading. 

At energies above the thermal peak, a power component can arise in the usual manner 
from shocks and synchrotron/IC radiation above the photosphere, and/or a comptonized 
tail in the photosphere itself. This approximation of the spectral peak determined by a 
black-body temperature _E p k ~ 3knT is only slightly changed if one considers the Compton 
equilibrium in dissipating photosphere which is optically thick to scattering, E p ^ ~ 3/cbT(1-|- 
A^ 1 ), where A ~ 1 is the photon to electron energy density ratio (Pe'er et al. 2005). 



2.2. Adiabatic Cooling and Radiative Acceleration of the Outflow 

This fireball radiation may be cooled adiabatically if it is generated close enough to 
the engine, so that the fireball Lorentz factor saturates at a value rj = Lj/Mjc 2 before the 
photons and the electrons are decoupled from each other. A condition for the neglect of 
adiabatic cooling is obtained by comparing the saturation radius _R sat with the photospheric 
radius. For the moment, we do not specify how outflow is accelerated and take 

r(r) = (-0 T (R T >r> Ro). (15) 

The saturation radius is then 

r~J Ro ' (16) 

If the electron-ion photosphere sits outside -R sa t, then its position is determined by setting 

aT W 2 (i^i^) = L (17) 
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Here fi is the mean mass per scattering charge. It is useful to defining the compactness 

cttL 



■'ISO 



87rm e c 3 r 3 ) _R 

at the thermalization radius, and the 'reduced' compactness 

r m P , „ 



(18) 



o = — *o- (19) 
/i 



Combining eqs. [IB] and [TTj gives 

p / _ \ (l+3a)/(l+2a) 

-rtsat _ V \ T-l/(l+2a) 



'0 



(20) 



Adiabatic cooling can be neglected if i? sat > -R T , which corresponds to 

v > r/cooi = t /(1+3a) r , (2i) 



or equivalently to 

This expression simplifies to 



-ftp L iso a T 



Ro n 1 n 9 -E"iso 52 

— > 6 x 10 j 

T tj !77| 



>6xl0 9 ^= cm (23) 



in the case of a ballistically expanding outflow (a = 1). Comparing this expression with eq. 
(EJ), one sees that adiabatic cooling can be neglect beyond the radius where the thermal peak 
energy is established, as long as rj = Lj/Mjc 2 > 10 2 . 

There is a close correspondence between the condition that adiabatic cooling be absent, 
and the condition that the thermal photon flux be strong enough to push the entrained 
baryons to a high Lorentz factor the electron-ion photosphere (see Grimsrud & Wasserman 
1998). One requires a mechanism for creating very high entropies in order to reach a terminal 
Lorentz factor as high as ~ 50-100. The photon field is quasi-isotropic in a frame that moves 
with a Lorentz factor 

r 7 (r) = r(i? r ), (24) 

where 

T(R T ) = 77 coo i I (rj > 77 coo i) (25) 

V^cool/ 
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is the combined Lorentz factor of the matter and photons at the photosphere. At r > R T , 
the radiation field has an anisotropic component in the rest frame of the matter, and the 
resulting force keeps T ~ T 7 out to a radius where 

ar 2T (i^v) ~ L (26) 

The matter therefore reaches a maximum Lorentz factor 

T^ = l l '\R T )T{R T ) ) (27) 

which is 



^?cool 

When a = 1, one recovers 

r^^A.^o^rr^)"' 4 ^)" 1 ' 4 . (29) 

The terminal Lorentz factor is smaller if the thermalization and acceleration are offset from 
the engine by a distance Rq 3> 10 6 cm. When a < 1, the terminal Lorentz factor is generally 
smaller than (1291) . and saturates this bound only if rj is large enough that the photosphere 
has contracted to R T ~ Rq. Note also that the entrained magnetic field does not limit 
acceleration by the thermal photons as long as T max (B 2 /Snpc 2 ) 1 ^ (Thompson 2006). 
Here p is the baryon density as measured in the rest frame of the central engine. 

This estimate of the terminal Lorentz factor allows for pair creation in the outflow out- 
side the baryonic photosphere. The conversion of a modest fraction of the outflow energy to 
pairs would significant decrease the inertia per scattering charge fi, and therefore increase the 
terminal Lorentz factor. In the absence of pair creation, the compactness at the electron-ion 
scattering photosphere is i ~ (m p /m e )(f]/T), which is easily large enough to allow effective 
collisions between photons of an energy ~ m e c 2 in the bulk frame. Only a tiny fraction 
~ 3/£ of the outflow energy must be converted to such photons to increase the number 
of light charges in the outflow. One can expect the relativistic ejecta to reach a terminal 
Lorentz factor as high as ~ 500 due to this effect. 

2.3. Implications of and Models for the Amati and Firmani Relations 

One can consider two basic solutions to eq. (jSJ) for the thermalization radius Rq, as 
deduced from the Amati et al. relation. In the first, the acceleration of the outflow is 



^?cool 



(l-a)/4(l+2a) 



cool J 



(2? 
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delayed, so that To ~ l/^o ~ few at a radius Rq. In the second case, the outflow has already 
attained a substantial fraction of its terminal speed, T ~ 77, when thermalization occurs. 

Case 1. When acceleration occurs after thermalization, one finds 



Here we have rewritten eq. (JSj) by substituting E iso = 2Ej/9j. 

A characteristic value of the thermalization radius follows from the following considera- 
tions. A certain fraction of the jet energy is mixed with the material of the Wolf-Rayet star, 
thereby forming a cocoon structure. When the energy in this cocoon becomes comparable 
to the binding energy of the Wolf-Rayet star, the star will explode and the rate of mass 
accretion onto the engine (the central black hole) must slow dramatically. 

After the jet head has reached the edge of the Wolf-Rayet core, relativistic material will 
continue to flow through the created opening. The Lorentz factor of this material generally 
increases away from the engine, but at a much slower rate than the Bernoulli rate for adiabatic 
flow. (This effect is observed in the simulations of an MHD jet propagating through a stellar 
envelope by McKinney 2006.) Two relativistic components of the jet can be distinguished: a 
central core which moves with a Lorentz factor T > 1/6; and an annular region sandwiched 
between the core and cocoon which is susceptible to strong non-radial shear instabilities. 
(These instabilities can involve several possible processes: oblique shocks, turbulent shear 
viscosity, and magnetic reconnection, the details of which are not addressed here.) 

This last jet component moves with a Lorentz factor F ~ 1/9; it is much hotter than 
the core but entrains only a small mass in baryons. Its width will be equated to the distance 
that a signal can propagate in the bulk frame, over the radial flow time Ro/2Tqc. This 
gives = 1/2-^/3, where we have taken the signal speed to be the sound speed cj ' \f?> in 
an isotropic relativistic fluid. Substituting this expression and Ej ~ 10 51 ergs (the binding 
energy of the core of a Wolf-Rayet star of initial mass ~ 25 M & ; Woosley & Weaver 1995) 
into eq. ( 1301) gives 



R = 1 x 10 10 (r ^) 




cm. 



(30) 



Ro 



3 x 10 




cm. 



(31) 



The peak energy scales with other quantities as 



E P k oc 



pl/4 ,1/4 T-, 

Ej tj R { 




(32) 



(Thompson 2006). 
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This interpretation of the Amati et al. relation has several interesting consequences: 

1. The thermalization radius f[3"Tj) corresponds closely to the radius of the CO core of a 
Wolf-Rayet star if the total jet energy is set equal to the core binding energy. The 
binding energy just before the collapse is approximately 1 x 10 51 (Mzams/20 M & ) in 
a star with a zero-age main sequence mass of 20-30 M (Woosley & Weaver 1995). 
The CO core binding energy and radius depend relatively weakly on details of the 
evolution of the progenitor - e.g. binary interaction, and the loss of a hydrogen or 
helium envelope. The typical duration of a long GRB is comparable to the collapse 
time of the core, whereas the light travel time across the core is only ~ 0.1 s. 

2. The temperature of the thermal radiation that is seen outside the photosphere is also 
insensitive to the baryon loading, as long as r\ lies below a critical value where tt > 
(k-BT^/rrieC 2 )' 1 at the thermalization radius. The electron-ion photosphere also lies 
well beyond the thermalization radius, 

p r 2 / 3 / p \ ~!/3 

W = 80 %F ( TnHT — ) £ 1051 er S s s' 1 )' ( 33 ) 

as long as the matter loading is high enough that the assumption of complete thermal- 
ization is justified. 

3. The relation between temperature and isotropic energy is insensitive to continuing 
dissipation over some distance outside the thermalization radius, if the outflow Lorentz 
factor grows linearly with radius (a = 1) 

r<r) = r„(£). (34) 

The scattering depth through the outflow decreases rapidly with radius, tt oc r _1 7~ 2 oc 
r -3 , so that the surface t t = {k-QTl h /m e c 2 )~ l sits at a radius ~ 0.2 — 0.4 R T . In other 
words, dissipation can continue over a factor 20 — 30 in radius beyond Rq and still 
result in a nearly black body spectrum; but continuing dissipation between this point 
and the photosphere will result in a broader 'greybody' spectrum (Pe'er et al. 2006). 

4. The thermal photons effectively accelerate the outflow even outside the electron-ion 
photosphere. As a result, portions of the jet that acquire higher entropies inside the 
Wolf-Rayet envelope may reach higher terminal Lorentz factors. A jet may contain a 
cooler core that flows with a Lorentz factor Y > 1/2^/36 and is subject only to weak 
shear-driven instabilities. This feedback between heating and acceleration provides an 
explanation for why a thermal radiation field should generically be present in the parts 
of the outflow that produce the prompt GRB emission. 
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Case 2. Now let us examine the other case, where acceleration occurs before thermalization. 
The temperature is now more sensitive to continuing dissipation in the outflow, and the 
relation between E pk and E iso is softer than implied by either the Amati et al. or Ghirlanda 
et al. relations. Suppose that the outflow Lorentz factor has saturated at T ~ rj. The 
thermal peak energy that results from dissipating a fraction e^b of the outflow energy at an 
optical depth tt — (&b^ / Vc 2 )" 1 is easily found to be 

= 2.7 k B T> = 0.5 4 7 b 6 J^j- {£) V3 keV. (35) 

In this case, the peak energy depends strongly on the baryon loading. In some models of the 
engine, e.g. those in which the outflow is driven by a magnetic field threading the horizon 
of a black hole (Blandford & Znajek 1977; Meszaros & Rees 1997), it is not obvious why r\ 
should correlate strongly with L iso . 

Recently, Firmani et al. (2006) have pointed out that a rather tight correlation ex- 
ists between three observed quantities that are derived from measurements of the prompt 
emission, namely L iso , E p ^ and to.45, namely L iso oc E pk 2 t^^ 9 . Here to.45 is essentially a 
measure of the duration of the prompt emission above a certain level (originally used for 
measuring the variability or spikiness of the prompt emission, by Fenimore & Ramirez-Ruiz 
2000, and Reichart et al. 2001). Assuming that t 0A5 ~ tj, and that the energy output rate is 
L iso ~ E iso /tj, the Firmani relation can be written as E- lso oc E^ k t l J 2 . The assumptions of a 
characteristic jet energy and thermalization radius, which were previously used to motivate 
the Amati et al. relation, imply 

E iso cx E 2 pk t) /2 Rl /2 (36) 

(see eq. [32] )• This essentially reproduces the phenomenological Firmani et al. (2006) 
relation. 



2.4. Correlation of E pk with Jet Opening Angle: Ghirlanda Relation 

Ghirlanda et al. (2004) derived a relation connecting the observed spectral peak energy 
E pk to the collimation-corrected jet energy Ej. The observed quantities are E pk , E iso and 
the break time tj, of the afterglow light curve. They are related by 

E pk cx EllX 2 , (37) 
(Liang & Zhang 2006, henceforth LZ; also Nava et al. 2006). 

A narrower distribution of burst energies results from the jet collimation correction in 
Ghirlanda et al. (2004), but a spread of about one order of magnitude in energy remains. In 
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this section, we relax our assumption of a single burst energy, and consider the implications 
of the Ghirlanda et al. relation for the jet properties within the Wolf-Rayet envelope. 

From t b and E iso one can obtain the opening angle 9j of the jet at the time of the 
afterglow break. We assume that the ambient medium has a power-law density profile, 

Pext oc r~ k . (38) 

The radius at which the break occurs is 

r(t b ) oc [£ iso /r 2 (t6)] 1/(3 - fc) oc 4 /{3 " fc) , (39) 
where T(t b ) ~ l/9j and Ej ~ (9 2 /2)E iso were used. Substituting 

t b ~ r(t b )/cT 2 (t b ) (40) 
into the observed LZ relation E^E^o, h) [equation ( (|37|) ]. one finds that 9j cancels out and 



<x Et k)/i6 - 2k) oc I E f ioik = (41) 
3 \ Ej for k — 2. V 1 



For a uniform external medium, k = 0, this yields the original Ghirlanda et al. relation 
i? p k oc E^ 3 , whereas for a wind, k = 2 it yields i? p k oc (a form which, as pointed out by 
Nava et al. (2006), is Lorentz invariant, since both quantities depend in the same manner 
on the bulk Lorentz factor). 

Taking the above correlation at face value, one can work backward and deduce a relation 
between the total jet energy and the opening angle of the jet. This relation is not unique, in 
that it depends on the index k of the density profile in the medium outside the progenitor. 
We focus, as before, on the interaction of the jet with the core of the Wolf-Rayet star, and fix 
the thermalization radius at the core radius Rq. Combining relation (T4T!) with the blackbody 
relation E pk oc \T(R ) / R^iE^/tj) 1 ^ and assuming the jet to be heated by shear-driven 
instabilities (T(R ) ~ 9^), one finds 

E J oc(9 2 t 1 / 2 R y 2 ^/ ( - 5 -^ (42) 

This gives Ej oc #~ 12 / 5 for k = and Ej oc #~ 4 / 3 for k = 2. The same scaling would, of 
course, hold if we chose some other fixed radius for the formation of the thermal spectrum. 

The Amati et al. and Firmani et al. relations both involve quantities that refer only to 
the prompt 7-ray emission. The Ghirlanda et al. relation combines quantities referring to 
both the prompt emission and the afterglow, which introduces an additional layer of model 
dependence. For example, the jet opening angle that is derived from the jet break time is 
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assumed to be the same as the opening angle in the prompt emission phase. A narrow range 
of outflow energies that is assumed in the derivation of the Amati et al. relation in §2.3[ 
whereas jet break correction suggests a somewhat broader range of energies (a factor ~ 10 
for a wind medium). The true range of jet energies is as yet undetermined. 



2.5. Heavier Baryon Loading and Possible Relation to X-ray Flashes 

Heavy baryon loading has two effects on an outflow: it tends to reduce the temperature 
of the thermal photons emerging from the flow; and it can also lengthen the duration of the 
outflow and hence of the photon signal. The first effect is generally encountered at lower 
baryon loadings than the second. 

X-ray flashes and GRBs are observed to have a similar distribution of durations (Sakamoto 
et al. 2005). Although the X-ray flashes are fainter and so could be observed from a smaller 
cosmological redshift, their intrinsic durations cannot differ by more than a factor 2-3 from 
those of the GRBs. One thereby obtains a strong constraint on the the baryon loading in 
the outflows that emit X-ray flashes. 

Consider an outflow in which the Lorentz factor has saturated at a value r ~ rj, and 
most of the energy is in the kinetic energy of the baryons. The photon signal is lengthened 
considerably if i] exceeds the critical value where 

h - (43) 

namely 

The photosphere generally sits outside the saturation radius (fTBl) when rj = rj t .. Substituting 
eq. (jHJ) into eq. (1201 and taking T = r/R gives 

pl/15 ,1/5 / N -1/15 
Rr 1Q W^fM . (45) 



This means that dissipation must continue beyond the saturation radius if a significant 
thermal photon energy flux is to be observed. The peak energy resulting from the dissipation 
of a fraction £ bb = | of the energy flux at a scattering depth tt ~ (k^T' / rrieC 2 )" 1 is obtained 
by substituting eq. (144")) into eq. (f3~5l) . 

E; k = 2.7 (tj) k B T> = 50 E}i\ 2 tjl /2 keV. (46) 
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Note that the dependence of E pk on E iso in this expression is weaker than in a blackbody. 
We conclude that values of E l pk as low as ~ 50 keV are consistent with burst durations of 
~ 10 s. 

Lower thermal photon energy fluxes and temperatures are possible if the dissipation 
occurs deeper in the outflow, and the thermal photons are adiabatically cooled. A simple 
relation between E p ^ and Ei SO is not expected in this regime, due to the sensitivity of the 
output spectrum to the location of the dissipation and the level of baryon loading. It should 
be emphasized that the adiabatically softened thermal photons can still act as the dominant 
coolant for relativistic electrons farther out in the outflow - even if their energy density is 
lower than that of the ambient magnetic field - in some models of electron heating (e.g. 
Thompson 2006). In this case, the thermal photons must be upscattered in frequency to 
absorb a significant portion of the outflow energy, and the effects of adiabatic cooling will 
be reversed. 



2.6. Implications for Long, Dim, Hard Outliers to the Amati Relation 

As pointed out by Nakar & Piran (2005), a significant number of long, dim and hard 
BATSE bursts fall, for any value of their unknown redshift, above and to the left of the 
nominal E pk vs. E iso Amati et al. relation. A relevant point of the present thermal model 
(Thompson 2006) is that the predicted E pk -E iso relation must transition from the nominal 
E pk oc E\Iq law to a flatter black-body type law E pk oc E^J* at low values of £i so , when 
9j gets larger than a critical value (corresponding to T(R ) ~ l/2\/30 ~ 1, or 9 ~ 20 
degrees (from the Frail et al. (2001) relation L iso tj(8 2 /2) = 5 x 10 50 ergs which leads to 
r(R*) = 1/2^36 = 0.9 Lift 1 /*). In fact, the redshift-localized bursts appear to cut off, in 
the Nakar and Piran (2005) simulations, at about the point where one would expect to see 
the transition, at E iso ~ 10 52 ergs. 

However, there could be other factors complicating the E pk -Ei SO relation below 10 52 
ergs. There could be a mix of different types of events, such off-axis, dirty fireballs that are 
adiabatically cooled, and off-center fireballs in which the reheating of soft bremsstrahlung 
photons is limited by pair annihilation and the saturation of the temperature at 30-50 keV 
(SN 1998bw type events). Thus, one should be cautious on interpreting the E pk values of the 
couple of very soft HETE-II bursts that are frequently plotted on the Amati et al. curve. 



-17- 



2.7. Low-energy Bursts Similar to GRB980425/SN1998bw 

Long bursts associated with detected supernovae (GRB/SNe) are a sub-class whose 
conformity (or not) to the Amati et al. relation is complicated. Some of these events 
satisfy this relation (e.g. GRB060218/SN2006aj, Campana et al. 2006) while others do 
not (e.g. GRB980425/SN1998bw and GRB031203/SN20031w). The gamma-ray output of 
many GRB/SNe is intrinsically low, and there is evidence pointing towards the presence 
of a transrelativistic (r ~ 1) outflow, e.g. GRB060218/SN2006aj, GRB980425/SN1998bw, 
GRB031203/SN20031w. Other GRB/SNe which are not particularly dim, such as GRB030329/ 
SN2003dh) may also have a transrelativistic component. 

This component could be identified with the ejection of a thin, fast shell during the 
breakout of a shock across the Wolf-Rayet photosphere (Colgate 1974; Tan, Matzner, & 
McKee 2001), which can supply up to ~ 10 48 ergs in transrelativistic material; or, in more 
energetic events, with a jet cocoon (Ramirez-Ruiz et al. 2002; Pe'er, Meszaros, & Rees 
2006). A third possibility is a shell of Wolf-Rayet material that is entrained at the head 
of a relativistic jet (Waxman & Meszaros 2003). Such a 'breakout shell' is susceptible to 
fragmentation, and in the case of a wide jet (opening angle tens of degrees) with a relatively 
short duration (tj < 10 s), its rest energy can approach the total jet energy (Thompson 2006). 
In this case, the two components will intermix within ~ 10 11 — 10 12 cm from the engine, and 
may not develop a large Lorentz factor. (The breakout shell is lighter and accelerates more 
easily in more focused jets with E iso > 10 52 ergs.) 

The radius at which gamma-rays are emitted from mildly relativistic ejecta depends on 
the ejecta mass and the density of the progenitor wind. There are various possibilities: the 
emission radius could be identified with the deceleration radius of the relativistic shell (espe- 
cially if the ejecta shell is light and the wind is dense; Tan et al. 2001); with the photosphere 
of the wind (Wang et al. 2006); or with the photosphere of the ejecta themselves. In the first 
case, the emission occurs at a time t em ~ E e ^V w /2T^ Q M w c 3 = 20 E C j iS V w $ (r c j /2)~ 4 M~L 4 
seconds in the observer's frame. (Here V w is the Wolf-Rayet wind velocity, M w = M ii ,_ 4 x 
10 -4 M yr^ 1 , and E e j is the ejecta energy and r e j o is its initial Lorentz factor.) In the second 
case, the wind photospheric radius is R T =i = Y e aTM w /4irV w m p in the absence of pair cre- 
ation, and the emission time is t em = R T= i/2T^c = 20 Mu,_ 4 /r 2 jl4, 8 seconds (for a wind com- 
posed of alpha elements, Y e — \ electrons per baryon). In the third case, the ejecta themselves 
become transparent to scattering at t em ~ (Y e a T E C j/16ir r| rripC 4 ) 1 ^ 2 = 2 x 10 3 eI^T'^ 2 
seconds (again neglecting pair creation and assuming Y e — | within the ejecta). (Note that 
in this last case, the ejecta are too heavy to be decelerated significantly by the Wolf-Rayet 
wind before they become optically thin, if the ejecta energy is as large as ~ 10 50 -10 51 ergs.) 
These timescales can be compared with the durations of GRB 980425 (t em ~ 30 s) and GRB 
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060218 (tem ~ 3 x 10 3 s for the prompt thermal X-ray component). 



Pair creation could evidently be important in the emission of GRB 980425, since Epk ~ 
100 keV was above the threshold for thermal pair creation (~ 30 keV in the bulk frame). 
We consider the simplest case of pairs in chemical equilibrium with a Wien gas of photons 
(density n'^ and mean energy 3/cbT 1 '). The pair density is related to the photon density by 



n', + n'_ 



„2\-l 



7T\l/2 

2. 



-3/2 



exp 



m e c 
kWT 



(47) 



Setting tt ~ {k B T' /m e c ) 1 (see eq. |14] ) for the limiting radius of thermalization gives an 
implicit relation for the rest frame temperature, 



t\ 3/2 



m e c z 
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kWT' 
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where 



87iT%m e c 3 r 



6 X 10 2 E,; 



j48 



30 s 



(48) 



(49) 



is the compactness in the emission zone. When the temperature is lower than this critical 
value, it is not possible for the thermal pairs to upscatter soft keV photons that are advected 
with the ejecta from the Wolf-Rayet photosphere. The energy of the Wien peak cannot be 
pushed any lower by sharing the thermal energy amongst a larger number of photons. 

In the case of GRB 980425, eqs. flU and (09]) give k B T' ~ 50 keV and E pk = § r ej x 
?>k B T' = 200 T e j keV. This is about a factor of 3 too hard for r e j = 2, but it should be 
recalled that upscattering of soft keV photons freezes out at an optical depth r T ~ 10. After 
the heating rate slows down, the photons entrained by the optically thick pair cloud would 
cool adiabaticallyEI The temperature could drop by by as much as r T 2 ^ 3 = 0.2 if pairs did 
not annihilate and T remained constant; in fact some deceleration and annihilation will take 
place, which have opposing effects on the output temperature. A burst duration of ~ 30 



seconds is obtained if the initial ejecta Lorentz factor is T 



ejO 



2 M 



-1/5 
'-4 • 



The X-ray spectrum of GRB 060218 was much too soft for pairs to contribute signifi- 
cantly to the optical depth. (The compactness £' is inferred to be only ~ 10 in the emission 



2 The specific heat of the photons dominates that of the pairs at temperatures well below m e c 2 //cb- We 
are interested in the regime where the baryonic rest energy is comparable to that of the photons, so the 
baryons provide inertia. They do not, however, provide a significant optical depth if £' < (m p /m e )(r]/T). 
The scattering depth through the neutralizing electrons can be smaller than unity at the point where the 
pairs freeze out, if the pairs are heated continuously starting at a large optical depth. 
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zone if the ejecta are mildly relativistic.) The high-energy tail of the X-ray spectrum could 
be explained by the Fermi acceleration of soft keV photons at the forward shock as it passes 
through the photosphere of the Wolf-Rayet wind, if r e j ~ 1 — 2 (Wang et al. 2006). This 
has the interesting implication that the pre-burst mass loss rate is quite high, M ~ 10~ 2 M & 
yr _1 , in order to place the wind photosphere at a distance ct em ~ 10 14 cm from the engine 
(see also Dai, Zhang, & Liang 2006). 

It also turns out that the radius ct em in GRB 060218 is close to the scattering pho- 
tosphere of the ejecta themselves, if they are mildly relativistic and have a kinetic energy 
of ~ 10 51 ergs. Therefore a related possibility is that the keV photons are warmed up by 
turbulence within the ejecta (which would be excited by the differential motion of the rela- 
tivistic ejecta with respect to the fragments of the breakout shell; see §3.1.11) . It should be 
emphasized that the black body radius of ~ 10 12 cm for the keV thermal emission (Campana 
et al. 2006) is much smaller than the photospheric radius that is inferred from the burst 
duration. This is consistent with the calculation of photon creation in §2.11 

It is possible that all GRB/SNe combine a transrelativistic outflow with a relativistic 
jet, the relative strengths of the two components contributing in varying ratios to the non- 
thermal 7-ray emission. However, the Amati et al. relation clearly relates to the relativistic 
component, whose presence is inferred in the great majority of classical bursts; whereas in 
events where the trans-relativistic component dominates it appears that the Amati et al. 
relation is generally not satisfied. Our analysis in §2.31 and 12.41 applies to bursts where the 
relativistic jet component dominates. 

2.8. Implications for the E pk -E iso Relation of Short GRBs 

The short GRB population offers a nice test of the idea that the peak energy is fixed 
by thermalization inside the outflow photosphere. Although the two burst populations have 
different progenitors, the outflow that produces them must have a very high compactness in 
both cases, and may be driven by essentially the same mechanism (e.g. a MHD jet). The 
short GRBs have systematically higher peak energies than is implied by the Amati et al. 
relation (Donaghy et al. 2006), but so far only GRB 050709 has both a well defined redshift 
and measured spectral peak energy (Fox et al. 2005; Villasenor et al. 2005). 

It should be kept in mind that the black body temperature depends more directly on 
the outflow luminosity than the total burst energy. To obtain the relevant thermalization 
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radius, we can make use of the scaling (eq. 

r f 1/2 

R oc i^-. (50) 

The isotropic energy of GRB 050709 is ~ 1 x 10 50 ergs and its peak energy is ~ 100 keV, 
some 10 times higher than would be expected based on the value of alone. However, 
the burst's T 90 duration is ~ 0.07s, more than 100 times shorter than a typical long GRB. 
Combining equation (150]) with eq. (jSJ) one obtains a thermalization radius of Rq ~ 1 x 10 10 



cm, similar to what we deduced for a long GRB with Ei SO ~ 10 52 ergs s 1 and tj ~ 10 s. 

A much smaller thermalization radius is deduced for the giant flare of 27 December 2004 
from SGR 1806—20. Here the peak energy was ~ 500 keV and the isotropic luminosity was 
4 x 10 47 ergs s _1 (Hurley et al. 2005; Palmer et al. 2005). The thermalization radius is 
therefore inferred to be ~ 50 km, about the Alfven radius expected for that luminosity and 
for a dipolar magnetic field of 10 15 G. 



3. Implications for the Emission Mechanism 

3.1. Peak Width-Photon Frequency Correlation 

The light curves of most gamma-ray bursts contain multiple pulses. A pulse is typically 
narrower when observed in a higher-frequency waveband: one observes the scaling 

St oc £-°- 5 (51) 

(Fenimore et al. 1995). 

A pulse results from dissipation within some part of the outflow. Its duration is limited 
by particle cooling as well as by the differential propagation time of the radiation across the 
emitting volume. We consider both of these effects in turn. 

A characteristic pulse width resulting from dissipation at radius r and Lorentz factor 
r is St ~ r/2T 2 c. Multiple pulses of this width would naturally result from dissipation well 
inside the radius at which the reverse shock wave passes through^ the ejecta shell (e.g. Sari 



3 This separation between an inner dissipative zone where internal shocks (or reconnection events) occur, 
and an outer dissipative zone where the afterglow is generated, is an artifact of the assumption that the 
gamma-ray emitting jet contains one dominant component. In fact, it is plausible that the outflow contains 
two components, one of which is much denser and slower. The slow component may be derived from a thin 
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Sz Piran 1997). This broadening effect is geometrical, however, and the resulting pulse width 
is to a first approximation independent of frequency. 

3.1.1. Implications for Beaming 

Our interpretation of the E^-E^q relation implies that photons above the spectral peak 
of a GRB are thermal photons upscattered in energy. The existence of a direct mapping 
between the spectral peak energy of a GRB and the temperature of the seed thermal photons 
then implies that photons near the spectral peak must be emitted by material moving in 
the direction of the observer. The same conclusion applies to higher-energy photons. The 
frequency dependence of the pulse profile could also, in principle, depend on the orientation 
of the the emitting material with respect to the line of sight (with softer photons being 
preferentially emitted off axis). Our model suggests that such orientation effects are of 
secondary importance in explaining the correlation (I5ip . 

This observation has important consequences for the emission mechanism. The scaling 
f l5T]) is suggestive of the cooling of relativistic particles: the characteristic frequency of the 
synchrotron or inverse-Compton photons emitted by an electron of energy 7 e m e c 2 is pro- 
portional to 7g, and the cooling time is inversely proportional to j e . The pulse width is 
controlled by cooling (over some range of frequencies) only if the light travel time across the 
radiating plasma is shorter than the cooling time at each (observed) frequency. This implies 
that the size of the cooling region is L'/ct f < (£7 e ) _1 , where i > 1 is the compactness in the 
bulk frame. 

One can consider dividing up each causally connected patch in the outflow (of size 
ct' ~ r/2r) into cells of size ~ ct'/jg. A single causal patch then comprises ~ 7^ such cells. 
If the emission is nearly isotropic in each cell, then the emission from any of them will be 
detectable. On the other hand, if the emission is beamed into a solid angle ~ 1/t|) then 
only ~ 7 e cells are visible to any observer. In this case, the pulse duty cycle is given simply 
by the fraction of cells that experience strong dissipation. 

Strong beaming of this type is not expected if the cooling particles are accelerated at a 
shock. Two mechanisms for beaming have been suggested in the case where the outflow is 
magnetically dominated: bulk relativistic motion of the magnetofluid (due to, e.g., current- 



shell of Wolf-Rayet material that is entrained at the jet head (Waxman & Meszaros 2006). This breakout 
shell is susceptible to fragmentation, and the shell fragments provide an attractive trigger for dissipation and 
gamma-ray emission as they drift backward through the faster, relativistic jet material (Thompson 2006). 
This process is completed at about 10 14 cm from the engine for E\ so ~ 10 52 ergs and tj ~ 10 s. 
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driven instabilities; Lyutikov & Blandford 2003); and heating of the motion of the light 
charges parallel to the background magnetic field, due to Landau damping of high-frequency 
Alfvenic turbulence (Thompson & Blaes 1998; Thompson 2006). 

An additional form of beaming is expected if the outflow carries a slower and denser 
component. The fragmentation of the heavy component occurs over a range of angular 
scales as it is accelerated outward by the momentum flux of the lighter relativistic fluid. At 
a distance of ~ 10 14 cm from the engine, the differential Lorentz factor between the two 
components has been reduced to Ar ~ 2 (Thompson 2006). A fraction of the seed thermal 
radiation will side-scatter off the heavy component (which remains optically thick out to 
~ 10 14 cm from the engine), and can provide an enhanced coolant for relativistic particles 
in the relativistic fluid. (See Sikora, Begelman & Rees 1994 for related considerations in the 
context of Blazar jets.) The inverse Compton radiation so produced will be preferentially 
beamed with respect to the side-scattered thermal radiation by a modest Doppler factor 

~ Ar. 



3.1.2. Inverse Compton Cooling and Other Cooling Mechanisms 

A simple relation between cooling time and photon frequency is expected if each high- 
energy photon results from a single upscattering. The cooling time in the bulk frame of the 
outflow is t' cool = 3m e c/4:'y e aTU! y , and the observed cooling time t coo i is shorter by a factor 
1/2T. Setting 

E, = 7 e 2 £pk (52) 
and making use of the Amati et al. relation (eq. pQ), gives 

^0) = 0.8 d£*H y^Y' 2 , (53) 

- C/ iso52 v 7 

(Here E iso is the isotropic energy of the prompt emission.) Pulses of a width of ~ 1 s must 
be emitted within ~ 10 15 cm from the engine, if the outflow Lorentz factor is close to the 
limiting value of ~ 10 2 (eq. [29]). Requiring that the cooling time (1531) be shorter than 
r/2r 2 c implies 

£3/4 , E vl/2 

r < 2 x 10 15 =^ 7 cm. (54) 

rjj^-i Viookevy 

Pair creation in the GRB outflow just outside the baryonic photosphere has a strong 
influence on the pulse duration. The prefactor T^ 3 could be as small as ~ 10~ 2 if the outflow 
became heavily pair loaded while being accelerated by the anisotropic pressure of the fireball 
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photons (eq. [29]). Substituting the upper bound fl54|) on r back into equation fl53l) for 
^cooi(IC), one deduces 

7^3/4 / J? \ X /2 

( »°' (IC)<001 (T^(100W) S ' (55 > 

Broader pulses can, of course, result from dissipation in the outflow outside the radius fl54|) . 
but their duration is limited by causality and not by IC cooling. 

A more energetic particle is required to emit a synchrotron photon of a given energy, 
which means that the synchrotron cooling time is correspondingly shorter. The equivalent 
result for synchrotron emission may be obtained by setting B' 2 /8tt = {eb/s^)U' and E' = 
0.3^heB' /m e c. One finds 

t cool (synch) in - 3 4 /8 ( zf* 

Wic) - 2xl ° ^'^f 1 ' (56) 

Pulses as broad as ~ 1 s, with the frequency-dependent width observed in GRB light curves, 
cannot easily be explained by synchrotron emission. 

A frequency-dependent pulse width would also result from the damping of bulk rela- 
tivistic motion by Compton drag, but in that case the scaling with photon frequency would 
be different. Suppose that a small blob of speed A/3 [Lorentz factor Ar = (1 — A/? 2 ) -1 / 2 ] 
is created in the outflow (e.g. by relativistic reconnection) . As it slows down to a speed 
Af3 < A/3o, the size of the blob increases to ~ (1 — A(3)ct' d , where t' dvag is the drag time at 
Lorentz factor AI\ The duration of the emission is therefore St oc t'/ AT 3 . Since the energy 
of the inverse Compton (IC) photons is E 1 oc Ar 2 , one finds the scaling St oc E^ 3 ^ 2 , much 
stronger than is observed. 

The frequency dependence of the pulse width at high and low frequencies deserves special 
scrutiny, for two reasons. 



1. The size of the emitting cells in the outflow is more or less independent of the maximum 
Lorentz factor 7 e of the heated particles. This means that observations of gamma-ray 
pulses at energies well above ~ 1 MeV may provide a direct diagnostic of the cell size: 
above some characteristic value of E y , the observed pulse width should saturate at 
some minimum value ~ L'/cT. 

2. Photons detected below the spectral peak have not been upscattered by particles with 
relativistic energies in the bulk frame. The pulse width should, as a result, have a 
weaker frequency dependence below the spectral peak. 
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We can summarize our conclusions as follows. First, the frequency-dependence of the 
pulse width arises most naturally from IC cooling of random particle motion. Since the 
cooling timescale is much shorter than the flow time in the emitting frame, the dissipation 
producing each pulse must be localized at some radius. This implies that the motion of 
the heated particles is strongly anisotropic; otherwise the emission on shorter timescales (at 
higher frequencies) would be smeared out by the differential propagation delay associated 
with shell curvature. 



3.2. Implications for the Mechanism of Particle Heating 

The temperature of the seed thermal photons is independent of radius (in the absence 
of adiabatic cooling). The peak energy resulting from IC scattering of the thermal photons 
therefore depends most directly on the mechanism by which the particles are heated: it is 
sensitive to the distribution of particle energies. When the heated particles have a power-law 
energy distribution that is cut off from below at r ) m _i n m e c 2 , it is clear that will remain 
close to the thermal peak energy only if 7 min ~ 1. 

A value of 7 m i n close to unity is easily achieved if the particles are heated continuously, 
and the outflow is at least moderately optically thick to scattering. Stochastic acceleration 
by turbulence (Dermer, Miller, & Li 1996; Thompson 2006) or by electrostatic acceleration 
in reconnection layers (Romanova & Lovelace 1992; Larrabee, Lovelace & Romanova 2003) 
are possible heating mechanisms. This implies that the high-energy tail of the prompt GRB 
spectrum is generated near the scattering photosphere. 

The synchrotron self-Compton process (Ghisellini & Celotti 1999; Stern & Poutanen 
2004) requires 7 m i n ^> 1. If the particle heating is truly continuous, then this condition can 
only be satisfied outside the scattering photosphere. In that case, 7 2 lin ~ r^ 1 oc rr 2 ^/L iso . 
The self-absorption frequency is given by hu Ba ~ T7g a (i? / /i?Q ED )m e c 2 oc L^ 14 T 2 / 7 r~ 5 / 7 , 
where 7 sa is given by eq. (fTUj) . The peak of the Comptonized synchrotron spectrum is at 

E^~^ n hu B ^ V T^L7^r^. (57) 

The explicit dependence on radius nearly scales out here, but the exponent of L iso has the 
wrong sign. (In this model, it is not clear how T and Lj SO are related, since the seed thermal 
radiation field must be assumed to be absent.) 

The continuous heating of seed thermal photons by a second-order Compton process has 
also been considered (Thompson 1994; Giannios 2006); but this mechanism does not provide 
an obvious explanation for the correlation (jSTl) between pulse width and photon frequency, 
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or the lags of soft photons with respect to hard photons that are observed in some GRBs 
(especially those with broad pulses, e.g. Norris et al. 2005). 



3.2.1. Shock Acceleration and 7 m i n ~ 1 

Shock acceleration is a rapid process: particles reach an energy 7 e m e c 2 on a timescale 
comparable to the gyroperiod •y e m e c/eB. Only modest Lorentz factors (j e ~ 1-10) are 
needed to create an IC spectrum extending up to ~ 10 2 E pk . At such low energies, the 
particle distribution is modified by cooling only at a large distance downstream of the shock. 

Under what conditions can shock acceleration generate a particle distribution with 
7min ~ 1, with a significant fraction of the outflow energy deposited in the non-thermal 
particles? Three basic requirements must be satisfied: first, a large fraction of the particle 
inertia in the outflow must be in light charges (electrons and positrons); second, a significant 
fraction of the light charges which encounter a shock must undergo Fermi acceleration; and, 
third, the momentum distribution of the of accelerated particles must have an index close 
to —2 (or harder). Although the behavior of a pair-dominated, relativistic shock is not fully 
understood (see Hoshino et al. 1992 for particle-in-cell simulations of electron-positron-ion 
shocks), these considerations are largely independent of such details. 

Pair creation in the outflow occurs primarily through collisions between photons (Cavallo 
& Rees 1978; Baring & Harding 1997; Lithwick & Sari 2001). These pair-creating photons 
have a characteristic energy ~ m e c 2 in the bulk frame, and are a byproduct of the cooling 
of charges of an energy 7 e ~ (Y m e c 2 / E^) 1 / 2 . Cooling is faster than pair creation, and so 
the equilibrium pair density can be determined by balancing the rates of pair creation and 
annihilation. Essentially the entire energy of the injected non-thermal particles is converted 
to a non-thermal photon spectrum by rapid IC cooling. 

To simplify matters, we assume that this process occurs over some range of radius in the 
outflow, so that the pair density is already close to its equilibrium value. For a high-energy 
photon index (5 > 2, the positron creation rate is n' + ~ (0.1 o-^n'^imeC 2 ) (Svensson 1987), 
where n!^(m e c 2 ) is the photon density at the threshold energy E' = m e c 2 . Balancing this 
with the annihilation rate n' + = — (3/8)o"Tc(n' + ) 2 in a warm pair plasma gives 



in' + n_ )m P <? 



" \m e c 2 J \Ym e c 2 



2*6^(13-2)1-^) ~ e nth (P - 2) (58) 



for rest energy density in the created pairs. The total photon energy density can be 
divided into a non-thermal tail carrying a fraction e nth of the total, and a thermal peak 
that has been boosted in energy by the cooling of the remaining thermal pairs. Now 7 min is 
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minimized if a significant fraction of the pairs are converted to a non-thermal distribution 
at the shock. The net energy density deposited in non-thermal pairs with particle index p is 
then 

' Ii:in -m e c 2 (n'+n'_)^e nth U' (59) 



(P-2) 

Combining this expression with fl58l) and the relation p — 2 = 2(j3 — 2) for rapidly cooling 
particles gives 

/ p \-C9-2) 

= 2 (rvSO ' (60) 

For example, if /3 = 2.2, _E pk = 100 keV and T = 10 2 then 7 min = 7. Even a modest 
departure of the photon spectrum from a —2 index implies that 7 m i n is large enough to give 
a considerable IC boost to the energy of the seed thermal photons. 

Softer high-energy spectra are observed in many GRBs (e.g. Preece et al. 2000). These 
spectra could still be consistent with a hard IC spectrum in the emission zone, if the high- 
energy photon flux were degraded by pair creation at a larger radius (e.g., in a dense Wolf- 
Rayet wind: Thompson & Madau 2000; Meszaros, Ramirez-ruiz, & Rees 2001; Beloborodov 
2002). 



4. Conclusions 

Most models of GRB outflows assume that they contain a single dominant component: 
e.g. baryons with a variable Lorentz factor (in the case of internal shock models) or a non- 
radial magnetic field with a reversing sign (in the case of reconnection models). Our central 
argument in this paper is that a second component is essential for understanding the prompt 
emission of GRBs: blackbody radiation that is emitted where the GRB jet forces its way 
through the core of a Wolf-Rayet star. Its thermal peak (Doppler boosted by the outflow) 
is identified with E pk . The non thermal high energy part of the GRB emission arises from 
Comptonization of this radiation by relativistic electrons outside the effective photosphere. 

This model accounts naturally for the small scatter in the Amati et al. and Firmani et 
al. relations. It should be re-iterated that both these relations can be tied directly to the 
jet properties in the thermalization zone, and are reproduced without free parameters by a 
very simple model. (Indeed, the relative dependence of E pk on the jet energy and duration 
that was measured by Firmani et al. 2006 was anticipated in the analysis of this model 
in Thompson 2006.) If we assume that the radius of prompt thermalization is fixed at the 
progenitor core radius (or some other characteristic radius), then it is possible to reproduce 
the Ghirlanda et al. relation for a particular dependence of jet energy on opening angle, 
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but not one that is motivated by a simple argument. A key, unresolved question centers 
on how the relativistic jet material builds up in the core of a Wolf-Rayet star: how wide is 
the outflow from the central engine, and how tightly regulated is the total energy that is 
deposited in the relativistic fluid with respect to the core binding energy? 

This quasi-thermal model for the observed spectral correlations has interesting impli- 
cations for the limiting jet Lorentz factor, and the relation of the jet energy, opening angle 
and burst duration to the mass and radius of the stellar stellar progenitor. The observed 
relation between pulse width and photon frequency can be explained by Compton cooling, 
but one requires that the relativistic particles in the outflow have an energy distribution with 
a low-energy (trans-relativistic) cutoff, and that the IC emission is beamed. The relation 
between E p k and E iso becomes more complicated for bursts with isotropic energies less than 
~ 10 52 ergs: the spectrum can be harder than is implied by the Amati et al. relation if the 
jet is clean but covers a large solid angle, or if the ejecta mass and energy are low. 
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